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QUARKONIUM PRODUCTION AND DECAY 

Richard S. Galik 
Cornell University, Ithaca, NY 14853 



ABSTRACT 

With new data sets, recently completed analyses and renewed interest, there has 
been significant progress in addressing existing questions about quarkonia produc- 
tion and decay ... but also new questions and new confrontation with theory. Some 
highlights include the firm establishment of the r]'^ and of the "missing" ip' decays 
to 0~1~ final states, improved information on the nature of the iIj{3770), and the 
observation of a new, puzzling charmonium-like state at 3872 MeV. 
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1 Introduction 



1.1 Why study Quarkonia? 

Quarkonia, bound states of a quark q and its anti-quark g, are the QCD equivalents 
of positronium (bound e^e~) in QED. Just as we have learned much about QED 
from the spectra of positronium, so too we probe the strong force with qq bound 
states. Quarkonia form the simplest, most symmetric strongly interacting system 
with only two constituents (unlike baryons) and those being identical (unlike mesons 
with "open" flavor). 

The QCD potential is much richer than that of QED. In the simplest 
formulation, the so-called "Cornell" potential, we have: 

V{r) = ---^ + k-r (1) 
6 r 

The first term mimics the single-photon exchange term in QED and corresponds 
to one-gluon exchange in QCD, dominating at short distances (< O.lfm) with large 
momentum transfers. This is the asymptotically "free" regime which has small, 
making calculations perturbative (PQCD). The second term, important at larger 
distance scales (> 1 fm), leads to "confinement" and is in a regime where is 
large, making calculations non-perturbative. 

In this presentation I will discuss only heavy quarkonia, QQ., with Q = c 
(charmonium) or h (bottomonium).^ 

A big advantage of such QQ systems is that they are not very relativistic: 
/32(cc) ~ 0.25 and P'^ihh) ~ 0.08. Therefore non-relativistic QCD (NRQCD) can be 
used to good approximation, greatly simplifying calculations. Of particular note^ 
is the fact that lattice QCD (LQCD) is making important advances in predictive 
power for these systems. Further, the physical sizes of QQ states span the range of 
"free" and "confined" , allowing both aspects of QCD to be studied. 

While this presentation concentrates on spectra and decays, one cannot 
ignore the importance of production in that much is to be learned from the formation 
of such QQ systems from (virtual) photons and gluons. 

The spectra for charmonia and bottomonia are shown in Figure ^ with 
some typical transitions indicated. Note that the QCD wells for such states are 
much deeper than the corresponding case for positronium, leading to a much richer 
spectrum, involving many hadronic transitions in addition to those associated with 

^Toponia would be wonderful to study, but the with mt > mw the top quark is not sufficiently 
stable to form hadrons before decaying. 

^See the contribution of Matt Wingate to this conference 
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photon emission. In particular, the bb spectrum has three radial excitations below 
open bottom threshold and enough energy splitting to allow multi-pion or u emission 
in transitions among states. 
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Figure 1: On the left (right) is the charmonium (hottomonium) spectrum below 
and near open flavor threshold. Shown for cc are di-lepton and two-photon decays, 
which, by t-reversal, are typical production mechanisms. The singlet h^. remains 
unconfirmed. Shown for bb are some of the transitions and decays with T ( 3S) as 
parent, indicating the diversity of modes and the large amount of available energy 
release, Q. None of the five singlet rjb or hb states have yet to be observed. 



1.2 Disclaimers and Previous Reviews 

This document is a faithful replication of the PIC04 presentation being a non- 
exhaustive selection of topics. Discoveries or updates made after PIC04 will not 
be included; this is therefore a snapshot of the field on 27 June 2004. One exception 
is that in the references I give the current status of the work. 

Of course this review builds on previous efforts, of which two good ex- 
amples are the PIC03 contribution of Mahlke-Kriiger and that from the 2003 
Lepton-Photon conference by Skwarnickij2]. There is also a comprehensive review 
of quarkonia underway by the Quarkonium Working Group this CERN Yellow 
Report is to be completed by September 2004. 

2 News on the QQ Spin-Singlets 

These states have the spins of Q and Q in an anti-symmetric state, so that S = 
and J = L. The states with J^*^' = 0^+ (pseudoscalars) are called "?7q" and those 
with J^^ = 1+- are called "/iq". 
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For bb there is no news on these singlets; the CLEO hmitsj2] from 2003 
have been shghtly updated[3] but none of the five have been observed as yet. In the 
cc sector these have been some new measurements of the mass, width and decays of 
the ground state of charmonium, the 7]^, but no starthng developments. The lone 
singlet P state of charmonium, the he, remains elusive, with active programs in pp 
production, ip' decay and B —>■ {cs)K (and maybe others!) all in the hunt. 

The biggest news from the spin singlets is the (re)discovery of the t]'^. This 
had been previously reported by the Crystal Ball experiment |4j in the Ml radiative 
transition ip' — > 777^. But CLEO, with similar sensitivity to the Crystal Ball,^ does 
not confirm that observation. [21 E] 

But four new observations of this radially excited singlet have now been 
published. First, Belle presented a very clean signal (see Fig. |2I) in B ^ rj'^K decays 
in]; they also observed the state in continuum production of double charmonium, 
e~^e~ J/ip r]'^ 17\. This was followed by both CLEOfS] and BaBar[n] publishing 
their measurements of the rj'^ in two-photon fusion. 

The four new measurements of the mass of the 1]'^ are shown in the right 
portion of Fig. |2l presented |10j as the hyperfine mass splitting between it and the 
spin-triplet ip'- The weighted average for this splitting is 

A'^^ = m(^') - m{r]'^) = (49 ± 2)MeV (2) 

which is roughly half the magnitude for this splitting based on the older Crystal Ball 
measurement of rn{ri'^, and which is to be compared with A/^j = m{J/ip) — mijic) = 
(117 ±2) MeV. ^ 

The splittings A/^j and A'^^ are of some interest because, while the rj^ and 
J /if) have small radius and are rather deep in the Coulomb-like QCD well, the r]'^ and 
start to sample the confinement region of the QCD potential. Most theoretical 
estimates of A'^^^ are higher than the new measurements [TU] . perhaps because they 
tend to assume a scalar QCD potential in the confinement region. A recent lattice 
result ■ which is quenched and thus ignores dynamical light quarks, gives a span 
of values of 40 < A'^^ < 74 MeV, depending on input parameters. This situation 
warrants a new look at models and an unquenched lattice calculation. 

^For example, CLEO does see the hindered transition to ^rjc at Str significance. 
''The four results have sufficient spread that using the PDGjJ prescription for scaling uncer- 
tainties leads to 49 ± 4 MeV for this hyperfine splitting. 
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Figure 2: First observations of the rj'^. On the left is the Belle spectrum in B ^ 
K{KgK7[) , showing the well-known rjc, and observation of the r]'^ at 3654 MeV. To 
the right is a summary of the four measurements by Belle, BaBar and CLEO. The 
weighted average of the hyperfine splitting to the ip' is 49 ±2 MeV. 



3 A Sampling of QQ Vector Results 

The vector ip and T resonances are directly produced by e~^e~ colliders via anni- 
hilation to a virtual photon. With large data sets recently obtained by CLEO at 
CESR and BES at BEPC, it is not surprising that these two experiments dominate 
this portion of the onia news. 

3.1 B^^ for the T States 

Good measurements of B(T{nS) = B^^ipS) are particularly important. 

They play a crucial role in determining the resonance widths, in that these states 
are significantly more narrow than the collider beam energy spread. CLEO plans 
to measure Fee to within a few percent from its careful scans of the resonant line 
shapes. Learning Ttot to ~ 5% thus means measuring to similar precision. Also, 
Bfj_^ is needed to evaluate many other important branching fractions, in that often 
what is really measured is the product B{X FT(nS)) ■ i3(T(nS)— > 

CLEO has both the statistical power (over 5 million of each of the three 
bound-state T resonances) and the control of systematics to make such precise 
measurements of B^^. Their preliminary results [T3] are i3^^(nS) = (2.53 ± 0.02 ± 
0.05)%, (2.11±0.03±0.05)%, and (2.44±0.07±0.05)% for n = 1, 2, 3, respectively. 
For the T(1S) this agrees with the current PDG averagejTTj, but the results for the 
two radial excitations are significantly higher than in the PDG. These results are 
shown graphically in Fig. El 
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Figure 3: Comparison of the new CLEO results for B^^ to prior measurements and, 
in the vertical stripe, the current PDC world average. 



3.2 From QQ to Q'Q' 

The color octet model was developed to help explain large charmonium production 
rates at hadron colliders by having cc produced from a single gluon (hence, as a color 
octet state) and then perturbatively shedding a soft gluon before emerging from the 
interaction. Models having two gluons form the cc meson as a final state particle 
are called color singlet models. A good place to test such octet [H] and singlet [T^ 
models is the glue-rich environment of T decay. While both predict roughly equal 
rates for T Jji\}X^ the singlet momentum spectrum is much softer in that a 
second c and c must hadronize as charmed mesons. 

CLEO HE] has measured i3(T ^ J/^X) = (6.4 ± 0.4 ± 0.6) x 10"^ which 
includes feed-down from other charmonia (■?/'', ...); this is consistent with the predic- 
tions of the models. However the measured J/?/' momentum spectrum is much too 
soft to support the color octet model. Further the observation of T — > ^/''X, with 
i3(T ^ i)'X) = (0.41 ± 0.11 ± 0.08) X B{T J/ipX), and the decay rates of T to 
Xci and Xc2 are above the octet model predictions. 

3.3 ^' V. J/i): The "12%" Puzzle 

Major progress has been made on a long standing puzzle in vector charmonium. Be- 
cause three-gluon decay to hadrons of cc and electro-magnetic production/decay of 
such states via a virtual photon both depend on the wave-function overlap (|^(0)p) 
one naively expects: 
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There are many complications, caveats and considerations ^ to this "equal- 
ity" , so one should not be too surprised at small deviations. However, two particular 
hadronic modes, pyr and K*K, both vector-pseudoscalar (V-P), are blatant violators 
of Eqn. El with limits [HI of Qh/Qee < 0.1. 

Big new data sets at BES/BEPC (14 million ip' decays and 6.4 pb~^ of 
continuum) and CLEO-c/CESR-c (~ 3 million ip' decays from 5.5 pb~^ of luminosity 
and ~ 20 pb^^ of continuum) have come to rescue. BES has presented evidence |17j 
for the modes K*'^K^ and K*^K~ in the final state KgK'^n^, as depicted in their 
Dalitz plot in Fig. |3J There is a large signal in the neutral mode, giving B{ip' 
K*°KO) = (15.0 ± 2.1 ± 1.7) X 10"^ and Qh = (3.6 ± 0.7)% (slight suppression 
relative to J/ip). The charged mode has a significant signal as well (3.5cr) and 
gives B{ij' K*+K~) = (2.9 ± 1.3 ± 0.4) x 10^^ and Qh = (0.6 ± 0.3)% (heavily 
suppressed relative to J/ip). Thus, there is a large difference in these isospin- related 
final states. The small continuum sample of BES does not allow for a subtraction 
of this possible background source. 
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Figure 4: On the left is the BES Dalitz plot from their analysis of ip' — ^ KgK^n^; 
there is a clear horizontal (vertical) enhancement corresponding to K*^Kg (K*^K'^) 
decays. On the right is the CLEO distribution of fractional visible energy in ip' — > 
TT+TT^TT^; the upper of the two spectra has had a p meson selected from the associated 
Dalitz plot projections. 



^For example, the running of as, form factor dependences on S, helicity conservation issues, 
non-relativistic effects, interference witli continuum, ... 
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CLEO has presented [TH] results for Qh/Qee for a large number of modes, 
including, for the first time, pvr in the tt+tt^tt^ final state. They have made a 
continuum subtraction, but assume no interference contribution. After projecting 
the Dalitz plot to obtain samples of p° and p^, they extract the yield by demanding 
no missing energy, as shown to the right of Fig. |31 Their results in the tt+tt^tt^ final 
state are Qp^/Qee = 0.016 ± 0.006 and Q^+^-no/Qee = 0.053 ± 0.011. 

3900704-012 

I I I 1 1 II 1 — I I 1 1 ii| r 

pjt 

(OJt°(IV) — 

pri(IV) I 

cor| — 

_ - 

K*OkO — 
K*+K~ 

cof2(1270) — 

^ cpfo(980) ^ 
fQ(980)— at+jt~ - 
cpf2(1525) — 

I I I L 

10~^ 1 10 10^ 10^ 

B(i|i(2S)) / B{J/^) (%) 

Figure 5: Qh, the ratio of branching fractions to specific hadronic states for the 
ip' = ip(2S) and J/ip- The solid vertical line is the naive expectation for this ratio: 
Qh = Qee = 12%. The inner dotted lines show the present uncertainties in Qee/ ihe 
outer dotted lines are at Qee/2 and 2Qee- 

The status of this "puzzle" is represented pictorially in Fig. El Some of the 
important features: (i) the pvr and K*K modes have now been measured] (ii) the 
isospin violating (and hence, electro-magnetic) modes seem to support the "12% 
rule"; (Hi) the isospin-related modes K*^K^ and K*^K~ seem quite different in 
their behavior; (iv) the A-P states do not seem suppressed, while V-T states are 
suppressed by about a factor of ~ 5 ... not nearly as much as the V-P states pvr and 
K*+K~. 

An open question is whether the suppression is due to S-D mixing, meaning 
a possible enhancement in V-P decays of the ^/^(STTO), by virtue of common virtual 
DD loops. 
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3.4 ip' V. J jil): Relative Phases 

Another interesting comparison of the Jjil) and wave-functions is the relative 
phase for each of these to decay strongly via cc annihilation to three gluons ("S") 
as opposed to electro-magnetically via a single virtual photon ("EM"). The decay 
(cc) K^K~ can proceed by either of these two routes, but (cc) K^K^ is purely 
a strong decay (SU(3) symmetry) and (cc) tt+tt^ is purely electro-magnetic (G- 
parity). Hence 

|^(ir+ir-)|2 = \A{^+^-)\''^A{K''W)\^+2-\A{^+^-)\-\A{K''W)\-cos<Ps,EM. (4) 

Previously BES[19J has measured cos(t)s,EM{,J l"^) = (90 ± 10)°. With 
14 X 10^ tp' events they have observed a large sample of mono-chromatic Kg can- 
didates from which they obtainjUI B^tp' K^s^l) = (5.24 ± 0.47 ± 0.48) x 10"^ 
Combining this with various prior results on the vr'^vr" and K^K^ channels gives 
cos(t)s,EM{i^') = (-89 ±29)° or (121 ±27)°. Thus, cos(Ps,EM{i^') and cos^s,EM{J/i^) 
are consistent. 

BES also re-measured |21j B{J/ip K^K^) from their 55 million J/ip 
sample as (1.82 ± 0.04 ± 0.13) x 10-^ meaning B{ip' K^Kl)/ B{J/ip KlKl) 
= (29 ± 4)%, somewhat large for the "12% rule" represented by Eqn. El 

4 QQ States with L = 2 ("D") 

Bottomonium (see Fig. is the only QCD system with a stable member that has 
two units of orbital angular momentum. This makes these "D" states important to 
test spin-orbit and high-L effects in both models and LQCD. The CLEO analysis 
of the four photon cascade T(3S)-^ 7172 "D" — > 71727374-^'''^" is now final P^. A 
single state is observed at Mo = (10161 ±0.6 ±1.6) MeV at > lOo" significance. The 
decay rate and intermediate x't and Xb assignments are consistent with this being 
the 1^D2 state. 

The impact of this result on LQCD|23j is shown in the left portion of Fig. [HI 
Plotted first results for nine "golden" quantities from analyses that were "quenched" , 
thus not incorporating the effects of dynamical light quarks. Next to this are recent 
itriquenched results with uds quarks included. None of the nine quantities are used 
as "input". These show LQCD being able to attain better than 5% accuracy in 
these quantities, one of which is the D-S splitting in bottomonium. 

Turning to charmonium, there is the continuing puzzle of the nature of the 
ip{S770), which one can characterize as being a ■ n^Si + (3 ■ l^Di ± 7 ■ DD. Clearly 
a 7^ in that the '?/'(3770) is directly produced in e"^e~ annihilation. Belle |24j has 
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Figure 6: On the left is the current status of nine "golden" quantities by which to 
compare lattice QCD calculations and experiment, with the most recent entry being 
the "D-S" mass splitting in bb. To the right is recoil mass against the two photons 
in T(3S)^ 77[7r+7r~T(^iSy]. The observed peak is due to El cascade to the T(2S); 
no rate is observed for decay from the new T{1^D2) state. 

recently observed a solid iIj{3770) signal in —>■ K^D^D'^, from which they extract 
B{B Kij{3770)) = (4.8 ±1.1 ±0.7) x 10"^ ~ 2/3 B{B Kip'). This comparison 
would imply that there is large S-D mixing, although color-octet models might be 
able to accommodate this rate for a pure "D" state (/? ~ 1). 

An important mode to investigate would seem to be ?/'(3770)^ -k^-k^ J /ip. 
A compilation of Markll and BES results|2S] gives T{'k+'k- J /i)) = (43 ± 14) keV, 
which is near the upper bound for this rate as analyzed by CLE0[2 . The Kuang-Yan 
prediction for this rate is 20-107 keV, depending on the level of S-D mixing; how- 
ever, the CLEO limit |22j . depicted in Fig. El for T(1D)^ ti^ti^J/iIj, is some seven 
times below 123 that predicted for a such a D state in the Kuang-Yan model, casting 
doubt on whether the T^ti^ti^ J/ip) measured for the '?/'(3770) is really consistent 
with theoretical expectations. Once the decay to -k^-k^ J/ip is solidly observed, the 
angular distributions of the decay products should help sort out|2H] 7- 

5 QQ States with L = 1 ("x") 

Although copiously produced in El transitions from the vector mesons directly ob- 
tained in e^e~ annihilation, rather little is known about these Xc and Xb states. 

CLEO has finalized j29j its analysis of x^j c<jT(1S) for J = 1,2. There 
is so little Q value for this decay that the process is kinematically forbidden for 
J=0! Nonetheless, the quoted branching fractions are large: B{x'bi ~^ ct;T(lS)) = 
(1.6±0.3±0.2)% and B{x'b2 ^ ^T(IS)) = (1.1±0.3±0.1)%, to be compared to the 
~ 7% branching fraction for the El photon transitions x'bj ~^ 7'^ (IS). The rates for 
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this oj transition are therefore roughly equal for J = 1 and J = 2, as predicted ^Hj. 

The BES collaboration has used its Xc sample to test the color-octet model 
in baryonic decays. The color singlet model is unable to generate the observed rates 
for, e.g., Xcj PP- Adding in color octet contributions jHO] predicts that the ratio 
TZb = r(AA)/r(pp) is ~ 0.60 and ~ 0.45 for Xci and Xc2 decays, respectively. The 
BES collaboration [31] last year reported on the AA channel; now they have very 
clean signals for all three J states decaying into pp. These lead to experimental 
values of TZb = 4.6 ± 2.3 and 5.1 ± 3.1 for J = 1 and J = 2. While the uncertainties 
are large, they nonetheless tend to show the AA channel enhanced, not suppressed, 
with respect to pp. 

6 The X(3872) - New Quarkonium-like State !! 

The hottest year's Lepton-Photon meeting(2] was the observation 

(> 10(t) by Belief! of a new "charmonium-like" state at 3872 MeV produced in 
K^X and decaying to J/ip ti^ti~ . Since then, it has been seen clearly by 
CDF|33j in hadro-production (as reported at QWGII^ and shown in Fig.[7j) as well 
as by DOPI and BaBar|S3]. 
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Figure 7: At QWGII, CDF confirmed the Belle observation of the X(3872), showing 
the left plot, with clear evidence of the ip' and the X (> llcr significance); the events 
were required to have a di-pion invariant mass in excess of 500 Me V. To the right 
are the BaBar spectra of this di-pion mass for both the X and ip' regions. 



Of particularly note is the nature of the tttt system in the X decay. Belief 
reported a dipion mass (m^^) consistent with the decay X p J/tp, which would 
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imply that Cx = +1- This tendency toward large m^^r "was actually used by CDF 
and DO to enhance their signals. Fig. [7| also shows the m^^r spectrum as reported 
by BaBar[35^ for the X and for their large '?/''(3686) sample. Is the high- mass region 
really a p or is it just a mimic of the similar shape in ip' decay? Is there a second, low- 
mass peak in the m,r7r spectrum, similar to what has been observed|3S] in T(3S)^ 
7r7rT(lS)? Do experiments see the more difficult channel X tt^tt^, which would 
imply Cx = — 1? One should keep watching this di-pion system as the sample sizes 
increase with more data at the Tevatron and the i?-factories. 

DOjni] has looked at six aspects of the production and decay of the X(3872) 
as compared to the '?/''(3686): to the jet, range of rapidity, isolation, decay length 
distribution, helicity angles of the muons and of the pions. In all cases the X looks 
like the well-established cc state. 

The four-experiment average has Mx = 3872.2 ± 0.5 MeV, tantalizingly 
close to the mass of a D^D^* pair at^T] 3871.2 ± 0.9 MeV but significantly below 
the corresponding charged system, with a mass of a D^D^* pair being 3879.3 ± 1.0 
MeV. This "coincidence", coupled with the fact [321 ES] that 

BiB+ ^ ir+X(3872)) ■ BjX ^ J/^ tt+tt-) ^ , ^ ^ 
B{B+ ^ i^+V^'(3686)) • Biij' ^ J/^ tt+tt-) ~ ^'^^ ^'^ ' 

leads to the conjecture that the X(3872) has a cc "core" and a large D^D^* "molec- 
ular" component. 

There has been a lot of activity (with a lot more ongoing!) to determine 
the quantum numbers of this new state. Searches by Belle [32j for the decays 7Xci, 
7Xc2, and 7J/?/' have shown it unlikely that X(3872) is l^Dz, V^D^, or Xc(2^Pj), 
respectively. Other studies of BES and CLEO data[37j in radiative return and two- 
photon fusion production tend to disfavor vector states ( J^*" = 1 ) and those with 
even J and positive C-parity. 

This is certainly a most curious state and hopefully much more will be 
known about it by PIC2005! 

7 Summary and Acknowledgments 

As is evident, heavy quarkonia continue to be a source of great energy and excite- 
ment, with many recent advances and many open questions: a stable "D" state 
has been firmly established and agrees nicely with the LQCD prediction; the long- 
awaited 7]'^ has also been firmly established with a small hyper-fine mass splitting to 
the several analyses are confronting the color octet production model; V-P states 
have finally been observed in ip' decay with suppressions relative to J/ip decay of 
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roughly a factor of 50; while more is being learned about the nature of the ■0(3770), 
it still remains a puzzle; even more of a mystery is the very narrow state at 3872 
MeV! Lots yet to do! 

I wish to thank all in the BES, BaBar, Belle, CDF, CLEO and DO groups 
who helped me put together this attempt at a selective summary. 
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